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ABSTRACT 

SN 1885 was a probable subluminous Type la supernova which occurred in the bulge of the Andromeda 
galaxy, M31, at a projected location 16" from the nucleus. Here we present and analyze Hubble Space 
Telescope images of the SN 1885 remnant seen in absorption against the M31 bulge via the resonance 
lines of Ca I, Ca II, Fe I, and Fe II. Viewed in Ca II H & K line absorption, the remnant appears as 
a nearly black circular spot with an outermost angular radius of 0'.'40 ± 0'.'025, implying a maximum 
linear radius of 1.52 ± 0.15 pc at M31's estimated distance of 785 ± 30 kpc and hence a 120 yr average 
expansion velocity of 12,400 ± 1400 kms -1 . The strongest Ca II absorption is organized in a broken ring 
structure with a radius of 0'.'2 (= 6000 km s _1 ) with several apparent absorption 'clumps' of an angular 
size around that of the image pixel scale of C/05 (= 1500 km s~ x ). Ca I and Fe I absorption structures 
appear similar except for a small Fe I absorption peak displaced O'.'l off-center of the Ca II structure by 
a projected velocity of about 3000 km s — 1 . 

Analyses of these images using off-center, delayed-detonation models suggest a low 56 Ni production 
similar to the subluminous SN la explosion of SN 1986G. The strongly lopsided images of of Ca I and Fe I 
can be understood as resulting from an aspherical chemical distribution, with the best agreement found 
using an off-center model viewed from an inclination of ~ 60° . The detection of small scale Ca II clumps 
is the first direct evidence for some instabilities and the existence of a deflagration phase in SNe la or, 
alternatively, mixing induced by radioactive decay of 56 Ni over time scales of seconds or days. However, 
the degree of mixing allowed by the observed images is much smaller than current 3D calculations for 
Rayleigh- Taylor dominated deflagration fronts. Moreover, the images require a central region of no or 
little Ca but iron group elements indicative for burning under sufficiently high densities for electron 
capture taking place, i.e., burning prior to a significant pre-expansion of the WD. Using time-dependent 
ionization calculations, we show that the presence today of neutral ions in this 120 yr old remnant can 
be understood as ejecta self-shielding from the UV radiation in the M31 bulge. 

Subject headings: supernovae: general - supernovae: individual (SN 1885) - ISM: kinematics and 
dynamics - ISM: abundances - supernova remnants 



The remnant of SN 1885 remained undetected for 
over a cent ury despite n umerous searches in t he optical 
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1. INTRODUCTION 

The 'bright nova' S Andromedae (S And) was discov- 
ered in the central bulge of M31 (NGC 224) during late 
August of 1885. It is believed to have been a Type la su- 
pernova (SN la) on th e basis of an apparent absence of hy - 
drogen Balmer lines (|Gaposchkinl ll936; Mi nkowski 1 19391 ) 
and a close match of nearly all reported spectral f eatures 
to those of SNe la (jde Vaucouleurs fe Corwinll 19851 ). How- 
ever relative to normal SN la events, SN 1885 was sublu- 
minous with an unusually rapid initial decline and a some- 
what redder B — V co l or (Ide Vaucouleurs fe Corwinl 1985t 
Patchett etal] 119851: iGrahaml 1 19881 IChevalier fe Plait! 
19881: Ivan den Berghlll994D. A recent analysis of SN 1885 
observations by Ivan den Berghl (|2002f ) suggests a decline 
rate of (1Mb (15) = 2.2 mag and a maximum brightness of 
M v = -18.74 and M B = -17.42, assuming a M31 dis- 
tance modulus of (m — M)o = 24.4 corresponding to a 
distance of 760 kpc. 

1 Based on observations with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope Science Institute, which is operated 
by the Association of Universities for Research in Astronomy, Inc. under NASA contract No. NAS5-26555. 

2 6127 Wilder Lab, Department of Physics & Astronomy, Dartmouth College, Hanover, NH 03755 

3 Department of Astronomy, University of Texas, Austin, TX 78712; and Department of Physics, Florida State University, Tallahassee, FL 
32306 

4 JILA and the Department of Astrophysical and Planetary Sciences, University of Colorado, Boulder, CO 80309 

5 Astrophysics Group, Imperial College London, Blackett Laboratory, Prince Consort Road, London SW7 2BZ 

6 Department of Astronomy and Astrophysics, University of Chicago, Chicago, IL 60637 
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The remnant was finally detected in 1989, not 
through its emission, but as a small dark spot of ab- 
sorption in the M31 bulge on an image taken using a 
near-UV filter (3900 ± 10 A) sensitive to C a II H & K 
resonance line absorp t ion j Fesen et"ail Il989l ) . Although 
ISiouwerman fe Dickell (|2001 ) report a possible weak radio 
detection at 8.4 GHz, SN 1885 continues to be essentially 
unobservable via emission, with recent deep Chandra X- 
ray imaging of the M3 1 bulge showi ng no associated X-ray 
source at its position (|Kaaretll2002l ). 

Following its ground-based detection, optical and UV 
imaging and spectra of the SN 1885 remnant (hereafter 
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A Ca II H & K absorption image of SNR 1885 was ob- 
tained with the ACS/WFC 2% ramp filter FR388N cen- 
tered at a wavelength 3950 A, chosen to match the center 
of the remnant's broad Ca II absorption line profile as 
observed i n the HST Faint Object Spectrograph (FOS) 
spectrum (jFesen et al.lll999T) . At this wavelength, the fil- 
ter's bandpass is 100 A FWHM with a peak transmission 
of about 50%. Images were taken on 2004 November 8 
with a dithering pattern and exposure times identical to 
those of the Ca I images. 

A continuum ACS/WFC image was also taken to allow 
for correcting small-scale variations in the distribution of 
background and foreground starlight from the M31 bulge 
at the position of SNR 1885. The continuum image was 
taken with the 9% ramp filter FR459M centered at a wave- 
length 4600 A, giving a bandpass of about 350 A FWHM 
and a peak transmission of nearly 80%. The spectral re- 
gion is adjacent to the Ca I and Ca II lines but free of any 
strong line absorption. Due to to the filter's wide bandpass 
and the increasing luminosity of M31 bulge stars to longer 
wavelengths, a relatively short exposure was adequate and 
dithered images were taken during just two orbits on 2004 
August 25 for a total exposure time of 4666 s. 

For obtaining UV Fe I and Fe II images, we used 
the available broadband filters on the higher resolution 
ACS/HRC. For the Fe Iimage, we selected the F330W fil- 
ter, which peaks at 3300 A and extends over 2975-3700 A. 
The filter catches t he Fe I reson ances lines at 3021 A, 
3441 A, and 3720 A (|Mortonlll99l[) . but most of its cover- 
age is line-free continuum. Consequently, the filter pro- 
vided only a weak detection of Fe I. Indeed, we used 
the Fe I image in part as a background continuum im- 
age to help correct the Ca I and Ca II images (§3SH1 E3- 
Dithered F330W images were taken over three orbits on 
2004 October 13 for a total exposure time of 7800 s. 

For the Fe II image, we selected the F250W filter, which 
peaks in transmission at 2500 A but has a broad trans- 
mission passband extending over 2250-3400 A. The filter 
encompasses the Fe II resonance lines at 2343 A, 2373 A, 
2382 A, 2586 A and 2599 A, but it also includes Fe I lines 
at 2524 A and 3021 A, the Mg I line at 2852 A, and the 
Mg II doublet at 2796 A and 2803 A, all o f which are 
predi cted to have significant optical depths (jFesen et alj 
1999). Dithered F250W images were taken over 6 orbits 
on 2004 November 20 for a total exposure of 15,600 s. 

Standard ACS pipeline IRAF/STSDAS^ data reduction 
was done, including debiasing, flat-fielding, geometric dis- 
tortion corrections, photometric calibrations, and cosmic 
ray and hot pixel removal. The STSDAS drizzle task was 
used to combine exposures in each filter. 

3. IMAGE RESULTS 

Figure [1] shows a 40" x 20" log scaled intensity view 
of SNR 1885 and its local environs within the bulge of 
M31. The Ca II H & K image (upper panel) reveals the 
remnant of SN 1885 as a dark and circular spot at the his- 
torically reported location of SN 1885. The coordinates for 
the center of the remnant are: a( J2000) = 00 h 42 m 42.96 s , 
<5(J2000) = +41°16'04"40. 

7 IRAF is distributed by the National Optical Astronomy Observatories, which is operated by the Association of Universities for Research 
in Astronomy, Inc. (AURA) under cooperative agreement with the National Science Foundation. The Space Telescope Science Data Analysis 
System (STSDAS) is distributed by the Space Telescope Science Institute. 



referred to as "SNR 1885") with instruments aboard the 
Hubble Space Telescope (HST) revealed a w 0"7 diameter 
absorption patch produced principally by Ca expa nding 
at velocities up to 13,000 km s _1 (jFesen et alJll999t ). The 
observed low flux at the Ca II H & K blended line center 
indicated a foreground starlight fraction of around 0.20, 
placing SNR 1885 some 65 pc to the near side of the M31 
bulge midpoint, comparable to its projected 55 pc (= 15"6) 
distance from the M31 nucleus. 

Other absorption lines detected in SNR 1885 included 
Ca I A4227, Fe I A3441, and Fe I A3720. A weak imaging 
detection of the remnant with the HST Wide-Field Plane- 
tary Camera-2 (WFPC2) in the wavelength range 2200 A 
to 3000 A, which is rich in strong Fe II lines, suggested a 
depth of absorption consistent with Fe II being fully satu- 
rated and an iron mass of 0.1 -1.0 M (jHamilton fe Fesenl 
I2000D . 

The fortuitous positioning of a relatively young and 
nearly freely expanding SN la remnant silhouetted against 
M31's bright central bulge offers an unique opportunity to 
map the distributions of Ca and Fe rich debris in a SN la. 
In this paper, we present images of SNR 1885 taken with 
the Advanced Camera for Surveys (ACS) on HST at a 
significantly higher angular resolution and signal-to-noise- 
ratio than previous images. The ACS images show a suffi- 
ciently well-resolved remnant to distinguish between vari- 
ous proposed SN la explosion scenarios. The observations, 
image results, our modeling results, future observations, 
and conclusions are given in §§2 — 6, respectively. 

2. OBSERVATIONS 

High-resolution optical and ultraviolet images of 
SNR 1885, detected via resonance line absorptions of M31 
bulge starlight, were obtained between 2004 August and 
November using the Advanced Camera for Surveys (ACS; 
iFord et alJH99alPavlovskv et al.ll2004h aboard HST. ACS 
has two main imaging configurations and both were used. 
The ACS Wide Field Channel (WFC) consists of two 
2048 x 4096 CCDs providing a field of view 202" x 202" 
with an average pixel size of 0"05, while the ACS High 
Resolution Channel (HRC) consists of a single 1024 x 1024 
CCD providing a spatial resolution of 0"028 x 0"025 per 
pixel and a nominal 29" x 26" field of view. 

The ACS/WFC was used in combination with narrow 
passband 'ramp' filters to image SNR 1885 in three sep- 
arate wavelength ranges. These were the resonance ab- 
sorption line of Ca I 4226.73 A, the Ca II H & K doublet 
at 3933.66 A and 3968.47 A, and a continuum band free 
of strong remnant absorption lines centered at 4600 A. 
Default ACS/WFC and ACS/HRC dither 'box' patterns 
were used which are 4-point parallelogram patterns with 
relative pixel coordinates (0, 0), (5.0, 1.5), (2.5, 4.5), (-2.5, 
3.0), thus having both integer and sub-pixel shifts. 

A Ca I image of SNR 1885 was taken with the 2% ramp 
filter FR423N centered at 4227 A, resulting in an effective 
bandpass of 75 A FWHM and a peak filter transmission 
of about 60%. Eight 1120 s exposures were taken during 
four HST orbits on 2004 November 2 in a four-point dither 
pattern for a total exposure time of 8960 s. 
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In contrast to the Ca II image, the offhand 4600 A 
FR459M continuum image (Fig. 1, lower panel) shows no 
appreciable absorption at the position of SNR 1885. Be- 
cause SN 1885's projected location is uncontaminated by 
the dust lanes that thread the M31 bulge, the lack of ap- 
preciable absorption seen in the 4600 A image at the po- 
sition of SNR 1885 suggests little internal dust formed in 
this SN la's metal-rich ejecta. 

Figure shows all four of the individual absorption im- 
ages of SNR 1885: Ca I A4227 (upper left), Ca II H & K 
AA3934,3968 (upper right), Fe I AA3441,3720 (lower left), 
and Fe II AA2358,2599 plus Mg II A2800 (lower right). 
The images shown here are in their native spatial pixel 
scales (ACS/WFC: 0'.'05; ACS/HRC: 0'.'028 x 0'/025) and 
are raw in the sense that they have not been corrected for 
any background variations in the surface brightness of the 
M31 bulge. 

The surface brightness of starlight from the bulge of M31 
has both a large scale gradient and smaller scale fluctua- 
tions. To correct the Ca I and Ca II absorption images 
of SNR 1885 for these non-uniformities, we used combi- 
nations of the offband 4600 A continuum image and the 
'Fe I' F330W filter image, in the ratio 4:1 and 3:1, re- 
spectively for the Ca I and Ca II images. As noted above 
in §2, the F330W filter is broad and, although it encom- 
passes a number of Fe I resonance lines, it mostly cov- 
ers unabsorbed continuum. We judged it better to use 
a background constructed from continuum both longward 
(4600 A) and shortward (F330W) of the Ca I and Ca II 
lines, and to accept the mild contamination from Fe I ab- 
sorption in the F330W continuum, than to use the 4600 A 
image alone. Because the ACS/HRC F330W image has a 
native image scale of 0'.'028 x 0'.'025 per pixel, we first re- 
sampled the ACS/WFC Ca I, Ca II, and 4600 A images to 
0'.'025 per pixel using the IRAF image interpolation task 
magnify, before applying the background correction. 

Figures [3] and [4] show in close-up the resampled, 
background-corrected images of SNR 1885 in Ca I and 
Ca II absorption, along with the original, uncorrected im- 
ages. We discuss these images in more detail below in 
§30 and O 

The distance to M31 is estimated to be between 760 
and 790 kpc based on studies of M31 Cepheid and 
RR Lyrae variables, globular clusters, and stars at the 
tip of the red giant branch (iFreedman fe Madorel 1990t 
Stanek fe Garnavich 19981: Joshi et al J 2003 : Brown et alj 
20041: iRich et al.ll2005t iMcConnachie et al.ll2005l and refer- 
ences therein). In the discussions below, we have adopted 
a distance of 785 ± 30 kpc to M31 and the SNR 1885. 

3.1. Cal 

Figure [3] shows the Ca I and 4600 A continuum images of 
SNR 1885 resampled to 0'.'025 per pixel. It also shows two 
slightly different background-corrected versions of the Ca I 
image. Like that seen in the raw Ca I image presented in 
Figure 2, the uncorrected but resampled Ca I absorption 
image has a markedly lopsided 'C shaped appearance with 
the strongest absorption oriented approximately in the di- 
rection of the nucleus of M31 toward the northeast. The 
remnant's overall diameter in this image is roughly 0'.'65 
corresponding to a linear diameter of 2.5 pc at a distance 
of 785 kpc. 



The two background-corrected Ca I images shown in 
Figure[3]are still noticeably lopsided, though much less so, 
and both appear more complete in the sense of showing 
a greater degree of absorption near the center and along 
the western limb. In the first background-corrected version 
(Fig. 3, lower left), we used only the 4600 A continuum im- 
age to remove bulge starlight variations, while for the sec- 
ond version (Fig. 3, lower right) which we regard as being 
better, we used a combination of 4600 A and F330W im- 
ages for the background in the ratio 4:1, respectively. One 
can see that the bulge non-uniformities are significantly 
decreased in both corrected images. Moreover, the good 
agreement between the two background-corrected Ca I im- 
ages indicates that the mild contamination from Fe I ab- 
sorption in the F330W image does not have a substantial 
effect. 

While our use of adjacent continuum images should 
largely remove small scale structures attributable to vari- 
ations in the background starlight, we cannot entirely re- 
move M31 bulge light variations because some 25% of the 
bulge starlight lies foreground to SNR 1885 (see M3.2I be- 
low). Nonetheless, the deepest levels of Ca I absorptions 
appear along the remnant's eastern and northern limbs. 

3.2. Call 

In the HST ACS images, SNR 1885 is most prominently 
visible through Ca II H & K line absorption. Figure 2] 
shows the ACS/WFC Ca II image both in its original 
0'.'05 pixel scale format and then after resampling to 0'.'025 
per pixel plus background corrected using the 4600 A and 
F330W images combined in the ratio 3:1. The intensity 
toward the center of the uncorrected Ca II image falls to 
about 30% of the background intensity of starlight, mak- 
ing SNR 1885 the darkest feature in the inner M31 bulge 
at 3950 A. 

Both original and background corrected images indicate 
the likely presence for a few Ca II absorption 'clumps' 
having an angular size around that of the ACS/WFC im- 
age resolution (0'.'05), ~ 0.20 pc at a distance of 785 kpc. 
This scale suggests a span in velocity for such clumps of 
around 1500 km s . The most prominent Ca II absorp- 
tion clump is about 25% darker than its surroundings and 
lies some O'/l due north of the remnant's Ca II absorption 
center. While some of the small-scale absorption struc- 
tures seen in the Ca II image are likely due to variations 
in both the corrected bulge starlight background and in the 
starlight lying in front of the remnant, the strength and 
non-randomness of the darkest Ca II features suggests the 
existence of a broken ring of Ca II-rich clumps at a ra- 
dius of 0'.'2 (6000kms _1 ) located along the inner edge of 
a limb-brighten shell. 

Aside from these apparent clumps, the degree of Ca II 
absorption in the image appears approximately constant 
from the center of the remnant out to a radius of 0'.'2, corre- 
sponding to a free expansion velocity of 6000 kms -1 . This 
suggests that the Ca II absorption is at least partially sat- 
urated in the central region, as previously con cluded from 
the FOS spectrum of SNR 1885 discussed by Fesen et al. 

con). 

The raw and background-corrected Ca II images of 
SNR 1885 also indicate a fairly spherical remnant struc- 
ture. The spherical shape and outer extent of the rem- 
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nant's Ca II absorption is illustrated by the 0'.'80 diameter 
circle shown in Figure 4. This circle encompasses virtually 
all of the observed absorption. 

The maximum detected dimensions of the SN 1885 rem- 
nant as seen in both the Ca II and Ca I images are shown 
in Figure [5l The Ca II image shows that Ca II absorp- 
tion extends out to a diameter of 0'.'80 ± 0'.'05, which is 
larger than that seen via Ca I absorption (diameter = 
0'.'70). Comparison of these Ca I and Ca II images which 
are shown here greatly stretched to bring out absorption 
in the outermost parts of the remnant, shows that the Ca I 
image as a whole appears offset relative to the center of 
the Ca II image, by approximately 0'.'05 in the direction 
of the nucleus of M31. The offset is consistent with the 
asymmetric appearance of the Ca I image, and probably 
reflects a real intrinsic asymmetry in the distribution of 
Ca I versus Ca II. 

The stretched Ca II image shown in Figure 5 indicates 
an outer radius for Ca II absorption of 1.52 ± 0.15pc at 
M31's estimated distance of 785 ± 30kpc. This is equiv- 
alent to an expansion velocity of 12,400 ± 1400 kms -1 
at the 120 yr age of the remnant. This velocity esti- 
mate is in agreement with the maximum Ca II velocity of 
13,100 ± 1500 km s" 1 directly measu red from a FOS spec- 
trum obtained byl^se^eFal. (1999). This agreement sup- 
ports the notion that the absorbing Ca II-rich ejecta are 
in free expansion. Moreover, since the image probes the 
size in the transverse direction while the spectrum probes 
the size in the line-of-sight direction, the similarity of the 
two expansion velocities indicates that the Ca II ejecta 
structure is approximately spherical. 

Finally, as noted above the average central intensity in 
the Ca II image is about 0.30 that of the local bulge inten- 
sity of starlight. This is somewhat higher than the aver- 
age central intensity of 0.25 seen in the FOS spectrum of 
SNR 1885 (|Fesen et al.lll999D . The slightly brighter aver- 
age central intensity measured from the ACS/WFC Ca II 
image compared to the FOS spectrum suggests that the 
foreground starlight fraction may be abou t 0.25 , higher 
than the 0.21 fraction found bv lFesen et al.l (1999) from a 
best-fit model of the FOS spectrum. A possible reason for 
this slight discrepancy in central intensities is that center 
of the 0'.'4 circular aperture of the FOS just happened to 
fall on one of the darker patches of Ca II absorption (see 
Fig. 4). 

3.3. Fel 

The remnant's Fe I absorption distribution was imaged 
using the broad, near-UV ACS/HRC filter F330W. While 
this filter's bandpass is sensitive to three Fe I lines (3021 A, 
3441 A, and 3720 A), the two strongest lines, at 3021 A 
and 3720 A, are detected only partially and at low filter 
transmission levels; that is, only the redshifted absorp- 
tion of the 3021 A line and just the blueshifted side of 
the 3720 A line. In addition, the F330W filter covered a 
considerable continuum region between 3000 — 3700 A out- 
side of the Fe I lines and, in fact, was used to help compose 
background bulge starlight images for the Ca I and Ca II 
ACS/WFC images. 

Nonetheless, faint Fe I absorption was detected on the 
F330W image. The Fe I image (Fig. 2) resembles most 
the uncorrected Ca I image, showing a lopsided and bro- 



ken absorption patch some 0'.'6 in diameter, a break along 
the southwest limb, and the strongest absorption concen- 
trated toward the northeast. The image also shows a peak 
of Fe I absorption roughly 2 — 3 pixels across (0'.'05 — 0'.'7), 
i.e., about the same size as that for the Ca II clumps, dis- 
placed toward the northeast by about O'.'l from the center 
of the Fe I absorption patch. 

3.4. Fell 

M31 bulge starlight is relatively faint between 2300 — 
3500 A and faintest between 2300 - 2600 A where several 
strong Fe II resonance lines are found. This fact, together 
with the rising continuum flux of bulge starlight toward 
the redward portion of the 2250 — 3400 A passband of 
the F250W filter on ACS/HRS resulted in a low signal-to- 
noise 3000 A continuum dominated image. This, in turn, 
produced a very weak detection of SNR 1885 presumably 
due to Fe II line absorption (Fig. 2, lower right panel) 
based upon ionization models which predict that the rem- 
nant should be almost black over the 2300 - 2600 A range 
(|Fesen et al.lll999t) . 

Notwithstanding its comparatively poor detection of 
SNR 1885, the F250W image bears a clear resemblance to 
the Ca II image of the remnant. This can be seen in Figure 
6 were we show a side-by-side comparison of SNR 1885 in 
Ca II and the F250W images. It is important to note that 
the F250W filter covers not only Fe II lines but also the 
Mg II doublet at 2800 A and several Fe I lines (see §3.3). 
However, the similarity of the F250W image to the Ca II 
image but not to the Fe I or Ca I images (see Fig. 2) sug- 
gests that Fe I does not make a major contribution to the 
F250W image. While Mg II 2800 A could also contribute 
to the detection of SNR 1885 on the F250W image, Mg II 
is expected to lie outside of the Ca-rich ejecta region and 
should thus lead to a larger diameter appearing remnant. 
Since this is not what is seen in the F250W image, we con- 
clude that the F250W image largely shows the distribution 
of Fe II-rich ejecta in SNR 1885. 

4. MODELS OF SNR 1885 

In this section, we compare our ACS imaging data of 
SNR 1885 to off-center delayed-detonation models of Type 
la SNe (SNe la). 

4.1. Overview of SN la scenarios 

There is general agreement that SNe la result from pro- 
cesses involving the combustion of a degenerate C-0 white 
dwarf (WD) that somehow reaches the Chandrasekhar 
mass limit, Mch- At present, two mechanisms are consid- 
ered most likely to bring a WD to this point. Currently, 
the leading idea is that a single WD accretes matter by 
Roche-lobe overflow from a main sequence or red giant 
companion ( Wh elan fe Ibenl Il973f ) . An alternative pos- 
sibility is the merger of two low-mass WDs in a binary 
(IWebbinkl fl98l llben fe Tutukov I [l98i iPaczvhskl 119851 : 
iBenz et all 119901 ) . Either way, the explosion is triggered 
near the center by compression. 

One of the main uncertainties in SN la physics is how the 
nuclear burning flame propagates through the white dwarf. 
Specifically, does it propagate as a weak supersonic deto- 
nation front, or as a subsonic deflagration wave? Whereas 
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detonation fronts are fairly well understood, modeling the 
propagation of a deflagration wave is difficult because the 
effective speed of deflagration burning depends on de- 
tails of the development of multi-dimensional instabilities 
which, in turn, depend sensitively on the initial conditions 
that lead to the thermonuclear runaway. 

Despite significant recent progress, direct deflagration 
modeling represents a major computational challenge. 
This is because of the multi-dimensional physics involved, 
the wide range of scales spanning 5 to 6 orders of magni- 
tude, and the sensitivity to initial conditions. As a result, 
fully self-consistent computer models are beyond reach at 
the present time, and additional observational constraints 
are needed to make progress on the subject. For example, 
observational knowledge of the distribution of the prod- 
ucts of nuclear burning, as in SNR 1885, can constrain the 
conditions under which burning must have taken place. 

Pure detonation explosions of Mch-niass models have 
long been ruled out for normal or subluminous SNe la. 
Such models predict that the WD is incinerated almost 
entirely to 56 Ni, contrary to observed SN la spectra which 
show a range of intermediate mass elements. Pure de- 
flagration explosions also appear to be ruled out for nor- 
mal SNe la. Three-dimensional numerical computations 
show that deflagrations produce deep Rayleigh- Taylor un- 
stable plumes, which lead to extensive radial , macroscopic 
mixing of burning products (lKhokhlovil200it iRopke et al.l 
12001 IGamezo et al.ll2005t iRopke fc Hillebrandtll2005D . As 



a consequence, pure deflagration models predict a chem- 
ically mixed region of burning products, surrounded by 
a thick layer of unburned C-0 of several tenths of a so- 
lar mass. By contrast, spectra of normal SNe la reveal 
a chemically stratified structure. Moreover, there is no 
evidence for rising deflagration plumes whose signatures 
should be visible in infrared spectra of normal SN la at 
late times o r in subluminous SNe la such as SN 199 9by at 
early times (jHoflich et al.ll2002t iHoflich et al.j|2004h . 

Better agreement with observations is obtained with so- 
called 'delayed-detonation' (DD) models in which the ex- 
plosion is assumed to begin as a deflagration which then 
(somehow) develops into a detonation. The mechanism 
by which a deflagration might turn into a detonation is 
not well understood. Thus simulations typically make the 
ad hoc assumption that at some point a deflagration wave 
simply turns into a detonation, and the density at which 
this transition takes place is treated as a free parameter. 

Increasing the transition density causes the transition 
from deflagration to detonation to occur earlier, increas- 
ing the amount of 56 Ni synthesized. Variation in this 
one parameter can account for the observed gross varia- 
tion in the properties of SNe la: how much 56 Ni is pro- 
duced, the overall distrib ution of elements, and the bright- 
ness of the SN la ev ent (|Hofjichlll995l : IHoflich et all 1 19951: 
llwamoto et~aLlll999l) . 

4.2. Delayed-Detonation Models 

Arguments about SNe la in general do not necessar- 
ily apply to an individual object such as SN 1885. How- 
ever, delayed-detonation models can reproduce both early 
to late time light curves and the spectra of both normal 
and subluminous SNe l a light spectra in the optical to the 
IR (IHoflich et alJl2003h . 



We therefore base our comparison to observations of 
SNR 188 5 on a suite of spherical delayed-detonation mod- 
els (jHoflich et al. 2002]). The models start from the pro- 
genitor, 5p01z22, a C-0 white dwarf taken from the 
core of an evolved 5 M main sequence star with solar 
metallicity. Through accretion, this core approaches the 
Chandrasekhar limit. A deflagration begins spontaneously 
when the core reaches a central densit y of 2.0 x 10 9 g cm" 3 
and the WD mass is close to 1.37 M© (jHoflich et al.ll2002D . 

As mentioned above concerning DD models, the density 
at which the deflagration is assumed to transition into a 
detonation is the key parameter that produces variation 
in SNe la. For extragalactic SNe la, the transition densi- 
tycan be adjusted so as to matc h the observed light curve, 
in particular its rate of decline (Hoflich et al.|[2003h . The 
chemical profile of the model then predicts a spectrum that 
can be compared to observation. In the case of SNR 1885, 
we will instead choose the transition density so as to repro- 
duce the chemical profiles observed with ACS. The maxi- 
mum brightness and light curve then become a prediction 
of the model. 

The HST ACS images of SNR 1885 presented in ig]show 
Ca II at velocities up to 12,400 kms" 1 . By comparison, 
DD models can produce Ca-rich ejecta with maximum ve- 
locities that range from 8000kms _1 for a very sublumi- 
nous SN la resembling SN 1991bg, up to 14,000 kms -1 
for normal SN la. The specific model that best matches 
SNR 1885 is 5p01z22.16, in which the '.16' refers to a DD 
transition density of 16 x 10 6 gcm~ 3 for the progenitor 
5p01z22. In this model, the deflagrated region encloses 
Mddt = 0.3 M©, and the explosion produces a net total 
of 0.27 M© of 56 Ni. 

Model 5p01z22.16, which also provides a good fit to 
the subluminous SN la object SN 1986G (jPhillips et al.l 
1987; Phillips 199^), is subluminous by about 0.5 magni- 
tudes, predicts a peak brightness of Mb = —18.72 and 
My = —18.74, and a 15-day decline rate of 5itlb = 1-70 
and 8mv = 1.26. This peak brightness in V agrees 
well with that observed for SN 1885, but the decline 
rate is slower than th e Sms = 2.2 reported for SN 1885 
(jvan den Berghl [20021 ). A possible explanation of the dis- 
crepancy in the decline rate may be related to uncertainties 
in the time t max of maximum light of SN 1885. All SNe la 
show a flat ma ximum, followed by a phase of increasing 
rate of decline (|Phillips et al.l[l9 99). A shift of —3 days in 
tmax would bring the model into closer agreement. 

4.3. The SN 1885 Explosion 

In constructing off-center D D models, we have followed 
the prescription of iLivnd (|l999f) . In this prescription, the 
initial deflagration phase is modeled assuming spherical 
symmetry. The deflagration begins at the center, and 
propagates outward in a subsonic deflagration wave. The 
energy deposited by deflagration causes the entire white 
dwarf to expand. When the density at the leading edge of 
the deflagration wave has fallen to a certain transition den- 
sity, a detonation is ignited by hand at a single point (the 
north pole) of the deflagration front. The simulation is 
then continued in 2 dimensions assuming cylindrical sym- 
metry. 

We choose to impose spherical symmetry on the initial 
deflagration phase for two reasons. The first reason is sim- 
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ply computational tractability: 3-dimensional simulations 
are co mputationally expensive, requiring many 10 5 CPU 
hours (|Khokhlovl 120011: iRopke et all 120031: iGamezo et all 
l2005h . and indeed no fully consistent DD models have 
ever been calculated. The second reason is that multi- 
dimensional simulations tend to predict mixed chemical 
profiles, which are at odds with observations both of typ- 
ical SN la (see JO above) and of SNR 1885, where the 
observed Ca I image indicates some radial stratification. 

No mixing seemed to us to offer a more promising ap- 
proach than full mixing. However, in H4.7I below we will 
show that the Ca image of SNR 1885 is more filled- in at 
its center than predicted by the model, suggesting that 
the assumption of spherical symmetry over-corrects the 
problem, producing more stratification than needed. We 
remedy this a posteriori in H4.71 by taking the final, fully- 
exploded WD, and partiall y remixing it s centr al regions. 

The 3D calculations of IGamezo et al.l (|2005[ ) show that 
detonation proceeds mainly along radial fingers, leaving 
behind connected regions of unburned matter. To allow for 
this, we assumed that the deflagration burns only about 
0.9 of the matter, leaving the remaining matter unburned. 
The unburned material plays an important role in that it 
allows the subsequent detonation to move back into and 
through the center of the partially deflagrated core, which 
offers a better appro ximation to the 3D calculations of 
IGamezo et al.l f|2005T ) than models in which there is no 
unburned material, and the detonatio n front is forced to 
circulate around the core iLivnel rfl999h . 

Deflagration burning causes the entire WD to expand. 
Several seconds into the deflagration, when the density has 
decreased by about two orders of magnitude, we initiate 
the detonation at a point on the boundary of the deflagra- 
tion wave. We treat the density of the transition as a free 
parameter, and we characterize its location by the mass 
Mddt interior to the detonation point. At the same time 
that we initiate the detonation, we switch from spherical 
to cylindrical symmetry, following the explosion on a grid 
of 256 x 90 cells in the radial and angular directions. The 
detonation begins at the pole of the cylindrically symmet- 
ric configuration, and propagates away from the pole in all 
directions. The detonation front moves not only outward 
and sideways through the unburned WD, but also inward 
through deflagrated core, consuming the 0.1 fraction of 
unburned material there. 

In the model, the detonation front reaches the surface of 
the WD at the pole about 0.2 seconds after the detonation, 
and it reaches the surface at the opposite pole about 0.4 
seconds after that, co nsistent with the tim es obtained in 
the 3D calculations of IGamezo et alj (|2005l ) (who did not, 
however, continue their calculation to homologous expan- 
sion). Over the next 5 to 10 seconds, the hot WD acceler- 
ates into free expansion, and the stored thermal energy is 
converted into kinetic energy. The density structure of the 
expanding ejecta freezes out as the sound speed plummets. 

About one second after entering the 2D phase, we switch 
back to ID hydrodynamics, evolving each of the 90 angu- 
lar segments as if it were spherically symmetric. At this 
time, angular gradients in the Ca-rich layers are less than 
2%. 

4.4. Numerical Methods of the Explosion Models 



All nu merical calcula tions were done using the HYDRA 
code (see iHoflichl l2003f) . HYDRA consists of a number 
of modules that solve simultaneously the hydrodynami- 
cal equations, nuclear and atomic networks, and radiation 
transport equations. 

For spherical geometry, the hydro equations are 
solved in the Lagrangian frame, and inc lude a 
front tracking scheme to re solve shock fronts ((H5flich 
I Wheeler fc Thielemann 1998f) . For cylindrical geometry, 
the coupled hydro and gravitational (Poisson) equations 
are solved in an Eulerian frame using a code based 
on PROMETHEUS (jFrvxellllArnett fc Muller 19911 ) with- 
out adaptive mesh refinement. The hydro modules 
use the explicit Piecewise P arabolic Method (PPM) of 
iCollella fc Woodward I (|1984l ) to solve the compressible re- 
active flow equations. 

HYDRA does not include a detailed model of the nu- 
clear flame and of the associated development of insta- 
bilities. Instead, the code parameterizes the deflagration 
speed semi-analytically, the parameters being adjusted so 
as to repr oduce deflagration speeds ob tained in 3D cal- 
culations. iDominguez fc Hoflichl (|2000l ) have shown that, 
in models with a parameterized deflagration speed, the 
abundance of chemical elements synthesized is relatively 
insensitive to the deflagration speed, but rather depends 
mainly on the amount of electron capture that takes place 
during deflagration, which in turn depends on the total 
amount of material burned during the deflagration phase. 

For the nuclear equation of state and the reactions, 
HYDRA us e s the nuclear reaction-netw ork library of 
I Thielemann! (|Nomoto fc Hashimoto 1994 and references 
therein). We take into account 218 isotopes for all mod- 
els. 

For spherical geometry, we use 912 radial layers, while 
for cylindrical geometry we adopt a relatively modest res- 
olution of 256 x 90 cells. Computationally, there is a trade- 
off between a large nuclear network and high spatial res- 
olution. We opt for a large nuclear network (rather than 
the 4 or 5 isotopes typically considered in 3D calculations) 
because numerical tests indicate that the final chemical 
composition is more sensitive to a reduction in the nuclear 
netw ork than to a factor 2 reduction in the spatial resolu- 
tion (jHoflichllWheeler fc Thielemann 19981) . In any case, 
the detonation phase is relatively brief and occurs on an 
already expanding WD so the effects of moderate spatial 
resolution during detonation are not expected to be dra- 
matic. Still, there is certainly room for improvement in 
this aspect of the modeling. 

4.5. The Chemical Distribution in the SNR 1885 Model 

Figure [7] shows the late-time chemical distribution in 
model 5p01z22. 16, with transition mass A/ddt = 0.3. This 
figure shows the distribution in velocity of O, Si, Ca, Fe, 
and 56 Ni several minutes after the explosion, which is well 
into the phase of free, homologous expansion. The det- 
onation axis points upward. Since the 56 Ni in SNR 1885 
today has long since decayed radioactively to 56 Fe, the fig- 
ure also shows the distribution of combined 56 Ni + Fe. The 
chemical structure seen is similar for models with transi- 
tion mass Mddt = 0.3 to 0.5 M Q . 

During the early stage of the deflagration, densities are 
high enough, > 10 9 gem -3 , that electron capture can take 
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place. This produces a central region containing stable 
isotopes of Fe, Co, and Ni with a mass fraction for Fe 
of 1=3 50%. As the density drops in the expanding WD, 
electron capture ceases and the later stages of deflagration 
yield about 0.1 M Q of radioactive Ni and Ca. 

The detonation phase produces off-center layers of ex- 
plosively synthesized isotopes, from 56 Ni and Ca on the 
inside to oxygen on the outside. The element abundances 
are characteristic of those obtained in various stages of 
burning. The outermost oxygen is partially burned, its 
mass fraction having increased from the 50% of the origi- 
nal C-0 WD to about 70%. 

The distribution of Ca in model 5p01z22.16 is qual- 
itatively similar to that seen in the HST/ ACS images 
of SNR 1885. The model shows Ca in a lopsided shell 
with a mass fraction around 3% from about 3500 to 
11,000 kms -1 , peaking at a mass fraction of about 6% at 
8200 km s -1 . The velocity of the peak increases slightly to 
8600 km s -1 if the transition mass is increased to A/ddt = 
0.5 M . 

The final Fe distribution has an inner region of about 
50%, produced during the early deflagration. This region 
is surrounded by a lopsided shell, produced by detonation 
to 56 Ni, where the mass fraction rises to almost 90%. The 
mass fraction decreases slowly outward, declining to about 
1% at 10,000 kms" 1 . 

The HST imaging observations of SNR 1885 reveal not 
the density itself, but rather the column density, i.e., 
the density projected along the line-of-sight. Thus Fig- 
ure [5] shows the column densities of Ca and Fe for Model 
5p01z22.16, with the detonation axis taken at various in- 
clinations to the line-of-sight. The projected distributions 
appear round when seen face-on (0°) but they become 
asymmetric at higher inclinations. The asymmetry re- 
mains modest even at high inclinations, in part because 
the intrinsic asymmetry is not large and in part because 
projection tends to smear out the asymmetry. Comparison 
between model and observed images of SNR 1885 does not 
tightly constrain the inclination. However, the observed 
Ca I image is clearly asymmetric, and we adopt an incli- 
nation angle of 60° because this angle is consistent with 
the projected asymmetry observed. 

4.6. The Ionization Structure of SNR 1885 

iFesen et al.l (|1999D and lHamilton fe Fesenl (|2000l) argued 
that the main source of ionization in SNR 1885 is photoion- 
izing UV from the bulge of M31, and that the main source 
of UV opacity in SN 1885 is from resonance lines of neu- 
tral and singly-ionized species. Recombination and charge 
exchange are negligible in SNR 1885 120 years after the 
e xplosion. 

IFesen et al.l |l999) also pointed out that if the ejecta 
were optically thin, then the lifetimes of Ca I and Fe I 
exposed to UV light from the bulge would be quite short; 
8-1 yr and I0~t\ 2 yr, respectively. The errors quoted come 
from uncertainty in the line-of-sight location of SN 1885 
within the bulge, but there could be a further factor of two 
uncertainty from photoionization cross-sections and from 
the extinction correction to the photoionizing UV flux. 

Predicted photoionization lifetimes as short as these 
pose a problem. Although resonance lines can block about 
half the photoionizing UV light at the present time, that 



still leaves half the UV to get through and red uce the Ca I 
and F e I below the currently observed levels. IFesen et al.l 
( 1999) suggested that Ca I and Fe I might be present, in 
spite of their short optically thin photoionizing times, if 
the remnant were optically thick in neutral continua, no- 
tably Fe I, earlier in its history. They did not, however, 
carry out the detailed model computations necessary to 
demonstrate this possibility. 

In the present paper, we have carried out the needed 
time-dependent, photoionization computations, and these 
show that the observed level of ionization in SNR 1885 
is entirely consistent with model expectations. We com- 
puted the photoionization of SNR 1885 as follows. We as- 
sumed an incident UV spectrum following the FOS data of 
IFesen et al.l (| 19991) do wn to 2221 A, from the IUE data of 
(|Burstein et"ailll988l) down to 1225 A, and from the HUT 
data of (jFerguson fe David son 1993) down to 912 A, dered- 
dened using the ICardelli et al.1 (Il989h extinction curve 
with color excess Eb-v = 0-11 (|Ferguson fe Davidseiil 
Il993h and R = A v /E B -v = 3.1. 

We treated the photoionizing UV spectrum incident on 
SNR 1885 as being isotropic. lHamilton fe Fesenl (|l99lh 
had noted that SNR 1885 should be receiving more radi- 
ation on the side facing the nucleus of M31 and suggested 
that for this reason SNR 1885 should be more ionized on 
that side. However, the ACS images show precisely the 
opposite. Not only do the images show Ca I concentrated 
toward the bulge, but an FOS spectrum of the remnant 
showed Ca I redshifted by 1 100 kms -1 , corres ponding to 
the side closer to the bulge (jFesen et al.lll999l ). Thus an 
anisotropic photoionization model cannot account for the 
observed asymmetry in the Ca I distribution. Rather, the 
observed asymmetry must reflect an intrinsic asymmetry 
in the explosion, such as that produced by an off-center 
delayed-detonation. For simplicity, therefore, we neglected 
any anisotropy of the incident photoionization. 

The evolution of the ionization structure of C, O, Mg, 
Si, S, Ca, Mn, Cr, Fe, and Ni was computed starting from 
about 10 years after the explosion. We took the initial 
ionization state to be mainly neutral, with an ionization 
fraction of 1% in the central deflagrated region where there 
is no radioactive heating, and 10% in the detonated region 
where there is radioactive heating, dominated initially by 

56 Ni _> 56 Cq _> 56p e and later by 44 Ti 44g c _, Uq & _ 

Although this initial degree of ionization is plausible, 
it is not well-determined a priori by the physics be- 
cause of complications associated with non-local energy 
transp ort by fast elec trons in the presence of magnetic 
fields dFransson et al.l Il996t iLiu. J efferv. fe Schultzl Il997t 



iRuiz-Lapuente fe Lucvl 1992 ). Tests indicated that the de- 
gree of ionization at the present time is most sensitive to 
the initial ionization fraction in the 56 Ni detonation layer. 
Reducing the ionization in this layer to 1% can delay ion- 
ization of Ca I and Fe I by about 30 years. The initial 
ionization state adopted here was chosen because it is rea- 
sonable and because it reproduces the observed degree of 
ionization at the present time. 

We computed the temperature structure by taking into 
account adiabatic expansion, radioactive heating, and the 
UV radiation field. The electron temperature affects re- 
combination rates and hence the ionization balance at 
early times when the density is high. 
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Line opacities were computed in the Sobolev approxima- 
tion for homologously expanding atmospheres, omitting 
stimulated emission. In homologously expanding ejecta, 
the radius equals velocity times age, r = vt, and the col- 
umn density per unit velocity dN/dv of any species is equal 
to its den sity n = dN/dr times the age t. In the nomen- 
clature of iMihalai (|l978l ). the Sobolev optical depth r at 
"wavelength A for a bound transition with absorption oscil- 
lator strength / is 

2 

r = fXnt (1) 

m e c 

where e and m e are the electron charge and mass, c the 
speed of light, and t is the 120 yr age of the remnant. 

Besides resonance lines, we also took into account ex- 
cited lines with lower levels originating from the ground 
level. We computed the level populations in a 2-level ap- 
proximation in which the lower levels are excited by (op- 
tically thin) IR and optical radiation from the bulge and 
which then decay by spontaneous emission. Altogether, we 
included 6,748 bound-bound transitions with wave lengths, 
energy levels, and oscillator strengths taken from iKuruc zl 
(|1995h . 

Both to follow photoionization and to take into account 
shielding by neutral atoms, we included bound-free opac- 
ities from the ground state for all neutral isotopes of the 
elem ents above. Photoion ization cross sections were taken 
from IVerner et al.l (|1996f ) to which we added some reso- 
nances treated in the narrow lin e approximation based o n 
data from the opacity project (jLvnas-Grav et all 1 19951) . 
Finally, we also took into account the Doppler shift of the 
continuum. 

The final result of all these calculations is that shielding 
of UV radiation by bound-free transitions of Mg I, Ca I, 
and Fe I proved to be the key in delaying photoionization 
so that Ca I and Fe I can survive in SNR 1885 up to the 
present time. All the models showed that Ca and Fe is 
mostly neutral at early times, but that Ca I and Fe I will 
rapidly disappear as soon as photoionizing UV radiation 
can penetrate the ejecta. Without bound-free shielding, 
Ca I is depleted to 10% after about 50 years, whereas with 
shielding, photoionizing radiation remains blocked for al- 
most a century. Whereas the UV opacity from bound-free 
transitions evolves substantially, most of the bound-bound 
transitions remain optically thick over the life of SNR 1885 
up to the present time, so that the UV line-blocking re- 
mains almost constant. 

Figure [9] shows the evolution of the Ca I/Ca fraction. 
The high-velocity, outermost layers having a solar-like 
abundance of Ca are quickly photoionized, whereas the 
inner, Ca-rich layers re main opaque for a bout a century. 
As previously argued bv lFesen et al.l (|l999l ). Ca-rich ejecta 
are currently undergoing a phase of rapid photoionization 
to Ca II. Indeed, our models predict that Ca I absorption 
will become unobservably small in just another 30 years 
or so, when the remnant is roughly 150 years old. 

4.7. 2D Model Images of SNR 1885 

Figure [TOl shows model Ca II, Ca I, and Fe I images that 
can be compared directly to the ACS images. These im- 
ages take into account not only the chemical and ionization 
structure of the remnant, as described above, but also the 
transmission curves of the HST filters, the contribution of 



different lines to absorption, and the finite angular resolu- 
tion of HST. As discussed in £|4.5i the detonation axis in 
the model is taken to be inclined at 60° to the line-of-sight. 

The model predicts the Ca II image to be insensitive 
to the precise ionization structure, since Ca is mostly 
singly-ionized. The model also predicts that the Ca im- 
age should have a sharp edge, where the mass fraction of 
Ca drops rapidly from 3% to the pre-explosion, solar-like 
abundance. The images shown in Figure [TU] are for model 
5p01z22.16, which as described in i )4.2l was chosen so as 
to reproduce the maximum velocity of Ca II observed in 
SNR 1885. 

The Ca II model image in Figure [TU] is as seen through 
the HST ramp-filter FR388N centered at 3950 land 100 
Awide, and takes into account the doublet character of 
the Ca II H & K absorption. The model predicts that the 
Ca II optical depth should exceed unity in a ring from 3100 
to 5500 km s -1 , which are close to the values measured for 
SNR 1885. 

The model Ca I and Fe I images in Figure [TUJ appear 
off-center by about 3000 and 2000 km s -1 , respectively, 
again in good qualitative agreement with the observed im- 
ages. The offset is caused primarily by the asymmetric 
chemical distribution produced by the assumed off-center 
delayed-detonation. The asymmetry is enhanced by the 
self-shielding of denser regions against photoionization. 

A defect of the model is that the predicted Ca I and 
Ca II images are emptier at their centers than observed 
in SNR 1885. It is natural to attribute this defect to our 
assumption that the early deflagration phase of the ex- 
plosion was spherically symmetric, which artificially sup- 
presses radial mixing by Ray lcigh- Taylor instabilities in 
the deflagration wave. Indeed, the observed images show 
structure down to scales ~ 1500 km s" 1 comparable to the 
angular resolution of the observations, such as might per- 
haps be produced by instabilities. 

To remedy this defect of the model, we tried radially 
re-mixing the final chemical profile in various ways. Fig- 
ure [TTJ shows the effect on the model Ca I and Ca II images 
of radially remixing of: (a) layers within the central de- 
flagrated region; (b) layers within the outer boundary of 
the Ca layer; or (c) the entire remnant. All three cases fill 
in the central Ca hole in the model to a greater or lesser 
degree. 

The observed Ca images of SNR 1885 are interme- 
diate between cases (a) and (b), suggesting that there 
may have been a moderate amount of radial mixing be- 
tween the inner deflagrated region and the surrounding 
detonated layers. This level of mixing is less than that 
predicted by curre nt 3D calculations of RT-dominated 
defla gration fronts (|Khokhlovll200H lR5pke fe Hillebrandt) 
2005). Smaller amounts of mixing could result if the ini- 
tial conditions at the ons et of runaway bur ning produce a 
less violent deflagration (|Livne et al.H2005h . or if the WD 
undergoes one or more puls ations with mixing prior to 
the transition to detonation ( 1 lollich et al.l l2~004). Mixing 
may also be produced or enhanced by heating of 56 Ni rich 
regions, well after the thermonuclear explosion, on time 
scales of hours to days. Both mechanism are compatible 
with fluctuations on scales observed. 

5. ALTERNATIVE SCENARIOS FOR SNR 1885 
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Because of the known diversity of SNe la, below we spec- 
ulate on some alternative scenarios for modeling the spe- 
cific case of SNR 1885. 

A variation of 'classical' DD models (which assume ad 
hoc that deflagration converts into a detonation), are the 
so- called 'confined-det onat ion' models r ecent ly suggested 
bv lCalder et"aH (|2003l) and lPlewa et~aH (|2004 . Their pro- 
posal is motivated in part by the extreme sensitivity to 
initial conditions of the early phase of burning. In this 
model, a buoyancy driven individual plume rises from the 
central deflagrating region to the surface, where it triggers 
a He-detonation. A detonation wave then propagates in- 
ward, exploding the WD. The confined-detonation model 
may offer a natural mechanism to burn a WD partly via 
deflagration and partly by detonation. 

Current models of confined-detonation use only a lim- 
ited nucleosynthetic network, so there is more work to be 
done, but the results to date suggest that burning in the 
central laye rs occurs at densit ies well below the value of 
10 9 gcm -3 (jPlewa et al.l 12004 ) needed to produce a cen- 
tral region with electron capture, such as is inferred from 
SNR 1885. Moreover, the confined detonation models pre- 
dict a significant amount of Fe at high velocities, which is 
not seen in normal SNe la. 

Mergers may contribute to the SNe la population and 
are obvious candidates for off-center explosions. For S- 
Andromeda, merger scenarios will suffer from the same 
problem as confined detonations namely, a low central 
density with ncgligable electron capture, unless the pri- 
mary star is close to Mch when the detonation is trig- 
gered. However, a detonation in a Mch mass WD would 
produce a large amount of Ni and consequently an over- 
luminous SNe la which is incompatible with observations 
of the SN 1885 event. 

Finally, one might even consider SN 1885 as a possi- 
ble core-collapse supernova. The fact that the remnant is 
freely expanding out to 13,000 km s _1 indicates that the 
remnant has so far swept up less than about 1.4 M© of 
material, implying an ambient circumstellar or interstellar 
density less than about 3 cm -3 . This relatively low den- 
sity, the historical record, and the absence of evidence of 
recent star formation in the bulge of M31, all point to a 
SNIa event, as opposed to a SNIb or SNIc event in which 
a massive star undergoes core collapse after having lost 
its hydrogen envelope in a wind. Our ACS observations 
add further evidence in this direction, the abundance and 
distribution of Fe and Ca being more consistent with a de- 
flagration/detonation event than with He-rich freeze-out. 
The near spherical symmetry of the Ca II image can also 
be construed as supporting the SN la interpretation, since 
there is growing evidence that the explosion mechanism of 
these objects is intrinsically asymmetric with axis ratios 
of 2: 1 or more (jWang et al.ll200ll iLeonard fe Filippenkol 
1200 lh . 

6. FUTURE OBSERVATIONS OF SNR 1885 

The remnant of SN 1885 offers an unique and remark- 
able opportunity to investigate the explosion dynamics 
and chemistry of a SN la explosion in some detail. The 
HST ACS data on SNR 1885 presented in this paper pro- 
vide only a glimpse of the potential to probe SN la physics. 
In this light, we thought it useful to consider future obser- 



vations that might provide further insights. 

Figure [12] shows a model UV absorption spectrum at 
three different times, respectively 90, 120, and 140 years 
after the explosion. The predicted spectrum is rich with 
strong, broad absorption features from dozens of resonance 
lines many of which evolve in strength fairly rapidly. 

Of course, the relative faintness of the bulge of M31 
in the UV will limit the observability of the SNR 1885 
spectrum. However, useful observations can be made even 
with limited signal-to-noise. For example, as discussed in 
N4.61 the remnant is undergoing a period of rapid ioniza- 
tion of Ca I and Fe I to the singly-ionized state, so that UV 
continuum absorption by neutral Ca and Fe is currently 
decreasing rapidly. This reduction in continuum opacity 
could be detected by broad-band UV imaging over a pe- 
riod as short as 5 years. 

Future spatially resolved UV spectral imaging, similar 
perhaps to that available with the Space Telescope Imag- 
ing Spectrograph (STIS), could provide a significant step 
forward in our understanding of SN la events. For ex- 
ample, Figure [T2] shows UV spectra taken with slits at 
various offsets from the remnant's center. Spectra such 
as these would provide a goldmine of information about 
the 3-dimensional structure and composition of the super- 
nova ejecta, and hence about uncertain aspects of explo- 
sion physics in SNe la. 

While narrow slits would of course provide the most 
precise mapping of ejecta, even broad slit spectra would 
provide useful information. Figure [T4l shows predicted UV 
spectra taken with various different slit widths available 
using the STIS on HST. 

7. CONCLUSIONS 

SN 1885, exploding as it did so close to, and to the near 
side of, the nuclear bulge of M31, presents a remarkable 
opportunity to map the chemical and density structure of 
freely-expanding ejecta in a Type la supernova. In this 
paper, we present high angular resolution Hubble Space 
Telescope images of the remnant of SN 1885 sensitive to 
Ca I, Ca II, Fe I, and Fe II resonance line absorption. 

The highest quality image is that of Ca II, which shows 
a spherical remnant 0'.'80±0'.'025 in diameter, correspond- 
ing to a free-expansion radius of 12,400 ± l,400kms _1 . 
Although the image is partially saturated (optical depth 
~ a few), it shows indications of a ring- like interior struc- 
ture peaking at a radius of about 6000kms _1 , with a few 
apparent absorption 'clumps' with an angular size near 
0'.'05 (1500 kms" 1 ). The SN 1885 remnant appears some- 
what lopsided and smaller (dia. = 0'.'65 ± 0'.'05) in Ca I 
than Ca II, indicating that the inner parts of the remnant 
are more neutral than the outer parts, as expected if the 
ejecta are being photoionized by UV light from the bulge 
of M31. The remnant's Ca I and Fe I absorption struc- 
tures appear similar in size and morphology, except that 
the Fe I image shows a small Fe I absorption peak of diam- 
eter 0'.'05 displaced to the northeast O'.'l (projected velocity 
of 3000 kms -1 ) from the center of the Ca II structure. Fi- 
nally, a weak detection of SNR 1885 across the wavelength 
range of 2300-3400 Ais probably caused by near saturated 
absorption from several resonance Fe II lines. The Fe II 
image appears similar in size to that seen in Ca II. 

We argue that the off-center and ring-like characteris- 
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tics of the Ca and Fe I images favor off-center delayed- 
detonation models of the explosion. We considered al- 
ternative scenarios including pure deflagration models, 
confined-detonations, and WD mergers, but find all of 
these to be less likely. 

Using the HYDRA code (Hoflich 2003), we have made 
detailed computations that follow delayed-detonation 
models of SN 1885 from the beginning of nuclear runaway 
up to the present 120 year age of the remnant. The mod- 
els follow in detail the nucleosynthesis in, radiative transfer 
through, and ionization structure of, the ejecta, including 
the effects of late-time photoionization by UV light from 
the bulge of M31. 

We selected a best-match model so as to reproduce the 
observed maximum velocity of Ca. The best-match model 
undergoes a transition from deflagration to detonation at 
a mass fraction of A/ddt = 0.3. The model is subluminous 
by about 0.5 magnitudes, consistent with indications from 
the historical record of SN 1885. The model synthesizes 
0.27 M Q of 56 Ni, and has a light curve similar to that of 
the subluminous Type la event, SN 1986G. 

We find that our model is able to reproduce many of the 
general observational features of SNR 1885. Ca and Fe in 
SNR 1885 are mostly singly-ionized at the present time, so 
the Ca II and Fe II images trace essentially all the Ca and 
Fe in the remnant. By contrast, Ca I and Fe I trace only 
the regions with the highest density and abundance of Ca 
and Fe, because only in these regions are Ca I and Fe I suf- 
ficiently self-shielding against photoionizing UV from the 
bulge of M31 to survive up to the present time. 

Our analysis presents the first direct comparison be- 
tween the properties of a supernova remnant and multi- 
dimensional explosion models to probe the underly- 
ing physics of burning fronts. Within the framework 
of delayed-detonation (DD) models, the deflagration- 
detonation transition occurs at densities of the or- 
der of 10 7 gm cm" 3 , consistent with the Zel'dovich 
mech anism (|Khokhlov et alj Il997t iNiemever fe Wooslevl 
|1997|) . Potential candidate processes for mixing of 
burned and unburned matter m ay be instabilities in- 
ternal to the deflagration front (Khokhl ov et al.l Il997t 



INiemever fe Wooslev 1997), during a pulsational phas e 
of the WD (KhokhloveialJ fl993 : H6flich_e±al| Il995h . 



or externally induces by shear instabilities present 
in rapidly, differentially rotating WDs when rising 
plumes enter this r e gion of strong differential rotation 
( Ho flich et all 120031: lUenishi, Nonioto, fe Hachisuj f2003: 
lYoon. Langer. fe Scheithauerll2004t lYoon fe Lang er 2004). 

The apparent presence of small scale Ca II-rich clumps 
is the first direct evidence for some instabilities and thus 
the existence of a deflagration phase in SNe la or, alter- 
natively, mixing induced by radioactive decay of 56 Ni over 
time scales of seconds or days. However, the degree of mix- 
ing allowed by the observed images is much smaller than 
current 3D calculations for Rayleigh- Taylor dominated de- 
flagration fronts. This could indicate a serious problem in 
the current generation of multi-dimensional models for de- 
flagration fronts, which may be related to the initial con- 
ditions at the time of the runaway. In addition, the images 
require a central region of no or little Ca but iron group 
elements indicative for burning under sufficiently high den- 
sities for electron capture taking place, i.e., burning prior 



to a significant pre-expansion of the WD. 

Off-center, delayed-detonation models have the virtue of 
being able to reproduce l ight curves and spe ctra of normal 
and subluminous SNe la (|Hoflich et al.|[2003T ). It is remark- 
able that these same models are able to reproduce the dis- 
tribution and morphology of Ca and Fe in the remnant of 
SN 1885. However, it should be emphasized that the mod- 
els described in this paper are parametrized, and thus not 
fully self-consistent. Our efforts should be regarded only 
as a first step toward tightening the link between explosion 
physics and observations of young supernova rem nants. 

Th e details of SNe la remain poorly understood fHoflich 
200(f) . How the deflagration begins and develops is ex- 
trem e ly sensitive to initial condition s (jHoflich fe Steinl 
I2002t iLivnej 12004 iGamezo et aTll2005l) . Equally unclear 
is how the deflagration subsequently turns into a detona- 
tion. While observations of any one object are unlikely 
to completely resolve such difficult questions, further ob- 
servations of the SN 1885 remnant might yield additional 
insights and guidance for future computer simulations of 
SN la explosions. 
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ported through UK PPARC grant PPA/G/S/2003/00040 
and JCW by NSF AST-0406740. 
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Fig. 1. — Upper panel: H8T ACS/WFC image of the M31 bulge centered 16 arcsec southwest of the nucleus. The ramp filter FR388N 
centered on 3950 A (FWHM = 100 A) was used to detect the remnant of SN 1885 (S And) via Ca II H and K absorption lines (3934, 3988 A). 
Lower panel: HST ACS/WFC of the same region but taken using the ramp filter FR459M centered at 4600 A (FWHM = 350 A). Comparison 
of these two images shows the presence of the SN 1885 remnant when imaged in Ca II. 




Fig. 2. — Upper panels: Drizzled ACS/WFC images of the SN 1885 in Ca I 4227 A and Ca II H & K at the native scale of 0'.'05 per pixel. 
Lower panels The remnant's appearance in Fe I and Fe II + Mg II line absorptions using the ACS/HRC (image scale: 0'.'028 X 0"025 per 
pixel) and the broadband filters F330W and F250W which are sensitive to Fe I and Fe II + Mg II, respectively. 
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Ofigiria! ACS/WFC Image 
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Fig. 3. — Resampled ACS/WFC images of SNR 1885 in Ca I (upper left), 4600 A continuum (upper right), a background corrected Ca I 
image resulting from a division of the Ca I image using the 4600 A image (lower left), and a Ca I image background corrected using a 4:1 
combination of the 4600 A and the F330W ACS/HRC images (lower right). See text for details. 



Original ACS/WFC Image 




Fig. 4. — Original ACS/WFC Ca II H & K image of SNR 1885 (left panel) and the resampled and background corrected Ca II (right panel) 
image The black circle shown is 0'.'80 in diameter. 
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Fig. 5. — Comparison of maximum Ca II and Ca I absorption diameters. Images are shown scaled as intensity squared with the contours 
of 0.85 and 0.97 for Ca II and Ca I, respectively. The Ca II and Ca I ACS/WFC images have been corrected for bulge background variations 
using the FR459M 4600 A (FWHM = 350 A) image and then resampled to 0'.'025 to match the higher pixel resolution of the ACS/HRC 
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Fig. 6. — Comparison of Ca II absorption with the Fe II ACS/HRC image taken with the F250W filter sensitive to Fe II absorption. 




Fig. 7. — Chemical distribution of 56 Ni and stable iron as a function of velocity for the off-center DDT at 0.3Mq. Boxi dimensions 
correspond to velocities from —27,500 to 27,500 km s -1 . The abundances are color-coded from zero (blue) to the maximum (red). The 
maximum mass fractions of O, Si, Ca, stable Fe and radioactive 5 ®Ni are 0.69, 0.65, 0.061, 0.59 and 0.90, respectively. In addition, we show 
the total Fe distribution about 120 years after the explosion. The central region consists of stable iron group elements Mn, Fe, Co, and Ni 
with mass fractions of ?s 4, 2, 58 and 36%, respectively (see Hoflich et al. 2002). The color bar is normalized to the maximum and used in all 
figures hereafter. 
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Fig. 8. — Column depths for the convolution of the density and chemical profiles for Ca and Fe as seen from various inclinations for the 
reference, off-center DD-model. Overall, the images are quite round. If viewed pole-on (90°), the chemical structure for Ca is a continuous 
ring which starts to break-up at lower latitudes. 




Fig. 9. — Cal/Ca as a function at year 40, 100 and 120 with a central value of 85%, 45%, and 16%, respectively. The ionization structure 
of Fe is rather similar in the models but, due to shielding by Ca I, larger by ss 3%. Note that there is no Ca I inside the 0.1 contour at day 
120. 
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Fig. 10. — Upper panels: Theoretical images Ca II, Ca I, and Fe I images of the off-center models seen at an inclination of 60°. The 
contours mark the layers with a darkness of 0, 0.25, 0.50 and 1.00. For Ca I, Ca II, and Fe I, the maximum darknesses are 0.64, 1.00 and 0.6, 
respectively. Lower Panels: Theoretical images convolved with the HST Ca II, Ca I, and Fe I filter response functions. Note that additional 
line blends tend to smear out the central hole due to the shift in frequency space. 




Fig. 11. — Influence of mixing on the Ca I and Ca II images in case of (a) mixing layers produced during the deflagration, (b) all layers 
within the Ca-shell, and (c) all layers that have undergone burning as expected for a (hypothetical) pure deflagration which burns the entire 
WD. For Ca I, case (b) and (c) are similar because Ca is fully ionized at the outer layers. 
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wavelength [A] 

Fig. 12. — Predicted UV spectra at 90, 120 and 140 years after the explosion with a spectrograph of a slit dimension of CK'l X 0'.'3 similar 
to that available on STIS. The slit has been assumed to be centered and aligned along the axis of symmetry defined by the off-center DDT 
and center of the white dwarf. Note the rapid spectral evolution below 2500 A. 





